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ABSTRACT 

We present new calibrations of far-ultraviolet (FUV) attenuation as derived 
from the total infrared to FUV luminosity ratio (IRX) and the FUV-NUV 
color. We find that the IRX-corrected FUV luminosities are tightly and lin- 
early correlated with the attenuation-corrected Ha luminosities (as measured 
from the Balmer decrement), with a rms scatter of ±0.09 dex. The ratios of 
these attenuation-corrected FUV to Ha luminosities are consistent with evolu- 
tionary synthesis model predictions, assuming a constant star formation rate 
over 100 Myr, solar metallicity and either a Salpeter or a Kroupa IMF with 
lower and upper mass limits of 0.1 and 100M o . The IRX-corrected FUV to 
Balmer-corrected Ha luminosity ratios do not show any trend with other galac- 
tic properties over the ranges covered by our sample objects. In contrast, FUV 
attenuation derived from the FUV-NUV color (UV spectral slope) show much 
larger random and systematic uncertainties. When compared to either Balmer- 
corrected Ha luminosities or IRX-corrected FUV luminosities the color-corrected 
FUV luminosities show ~ 2.5 times larger rms scatter, and systematic nonlinear 
deviations as functions of luminosity and other parameters. Linear combinations 
of 25/iin and 1.4GHz radio continuum luminosities with the observed FUV lu- 
minosities are also well correlated with the Balmer-corrected Ha luminosities. 
These results provide useful prescriptions for deriving attenuation-corrected star 
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formation rates of galaxies based on linear combinations of UV and IR or ra- 
dio luminosities, which are presented in convenient tabular form. Comparisons 
of our calibrations with attenuation corrections in the literature and with dust 
attenuation laws are also made. 

Subject headings: dust, extinction - - galaxies: ISM - - infrared: galaxies - 
ultraviolet: galaxies 

1. INTRODUCTION 

The ultraviolet (UV) continuum of galaxies traces young massive stars, and hence serves 
as one of the most commonly used star formation rate (SFR) indicators (e.g., Kennicutt 1998, 
hereafter K98). However, the severe dust attenuation in the UV significantly weakens its 
power for measuring the SFR. In nearby star-forming galaxies, the dust attenuation in the 
UV can vary from zero to several magnitudes (e.g., Buat et al. 2005; Burgarella et al. 2005). 
For near-UV (NUV) selected and far-infrared (FIR) selected samples, Buat et al. (2005) 
found that the median dust attenuations in the far-UV (FUV) are 1.1 mag and 2.9 mag 
respectively (see also Burgarella et al. 2005), corresponding to 64% and 93% obscuration of 
the FUV light respectively. Corrections for dust attenuation therefore are essential before 
the UV emission is used to probe the SFR. 

Commonly two approaches are used to estimate the dust attenuation in the UV. The 
power-law slope of the UV continuum (3 (fx oc X 13 , at wavelengths 1300-2600A; Calzetti 
et al. (1994)) was the first proposed UV dust attenuation indicator (e.g., Kinney et al. 
1993; Calzetti et al. 1994; Meurer et al. 1999), and is based on the assumption that the 
intrinsic slope /3 is only sensitive to the recent star formation history and any deviation 
from the intrinsic value is caused by dust attenuation (e.g., Calzetti et al. 1994; Meurer 
et al. 1995). However, once similar analysis is extended from starburst galaxies to normal 
star-forming galaxies, j3 will depend on the mean age of the dust-heating population, the 
dust/star geometry, and the dust properties (e.g., Witt, Thronson & Capuano 1992; Witt & 
Gordon 2000; Granato et al. 2000; Burgarella et al. 2005). 

The other UV dust attenuation estimator is the ratio of IR to UV emission. This 
method is built on an energetic budget consideration (e.g., Meurer et al. 1999; Buat et al. 
1999). According to the energy balance argument, all the starlight absorbed at UV and 
optical wavebands by the interstellar dust is re-emitted in the IR, so the combination of 
the IR luminosity and the observed UV luminosity should be able to probe the dust-free 
UV luminosity (e.g., Wang & Heckman 1996; Heckman et al. 1998), which then can be 
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used in turn to estimate the dust attenuation by comparing to the observed UV luminosity. 
Compared to (3, the IR/UV ratio is a more reliable dust attenuation indicator because it is 
almost independent of the dust properties and the relative distribution of dust and stars (e.g., 
Buat & Xu 1996; Gordon et al. 2000; Witt & Gordon 2000; Cortese et al. 2008). However, 
the IR/UV ratio does depend on the age of the dust heating populations (e.g., Cortese et al. 
2008; B. D. Johnson et al. 2011, in preparation). The non-ionizing UV emission is mainly 
from stars younger than a few hundred million years, whereas by contrast the IR emission 
can be from dust heated by stars over a wide age range, including not only young stars that 
emit a bulk of their bolometric luminosities in the UV but also evolved stars that radiate 
only a small fraction of their bolometric energies in the UV or even do not contribute to the 
UV emission at all. 

Since the introduction of these two methods many efforts have been made to calibrate 
the relations between UV attenuation and the IR/UV ratio or j3 (or equivalently FUV-NUV 
color) or the IR/UV versus /3 relation (e.g. Meurer et al. 1999; Buat et al. 1999; Gordon et 
al. 2000; Bell 2002; Kong et al. 2004; Calzetti et al. 2005; Seibert et al. 2005; Cortese et al. 
2006; Boissier et al. 2007; Gil de Paz et al. 2007; Johnson et al. 2007a,2007b; Panuzzo et 
al. 2007; Salim et al. 2007,2009; Treyer et al. 2007; Boquien et al. 2009; Buat et al. 2010; 
Takeuchi et al. 2010). Most of the calibrations are either purely based on theoretical models 
or based on comparisons of mult i- wavelength photometry with stellar population synthesis 
models. For example, Meurer et al. (1999) obtained a relationship between FIR/FUV and 
the attenuation at FUV (A FU v) using the stellar population synthesis models of Leitherer & 
Heckman (1995) and then they derived an A F uy-/3 relation using this FIR/FUV-based A FUV 
for starburst galaxies. Buat et al. (2005) derived Auv~IR/UV relations for star-forming 
galaxies at FUV and NUV wavebands based purely on stellar population synthesis modelling. 
Salim et al. (2007) estimated the UV attenuation by comparing the observed spectral energy 
distributions (SEDs) to modelled SEDs and then based on this UV attenuation they obtained 
A F uv~f3 relations for star-forming galaxies selected spectroscopically or using NUV-r color. 
Other studies that use different samples or types of galaxies to derive prescriptions for 
estimating UV attenuation either from IR/UV or from f3 (or equivalently FUV-NUV color) 
include Burgarella et al. 2005, Seibert et al. 2005 and Cortese et al. 2008. As discussed 
earlier, the IR/UV ratio and the FUV-NUV color are dependent on the age of the dust 
heating populations, so calibrations based on different samples, regardless of the methods 
adopted, should account for this age effect. The considerable inconsistencies between the 
above prescriptions (see Section 6.2) imply that the age difference of the stellar populations 
of the samples used was largely taken into account, and these prescriptions have been applied 
widely in the literature to correct for UV attenuations (e.g., Bell & Kennicutt 2001; Iglesias- 
Paramo et al. 2006). However, the reliability and accuracy of these corrections need to be 
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evaluated using an independent set of attenuation estimates. 

In the literature an independent attenuation estimate, the attenuation in the Ha line 
(A Ha ) as derived from the Balmer decrement ratio, has also been used to estimate the 
attenuation in the UV. To achieve this, two methods have been employed. One is based on 
the idea that the SFR derived from attenuation-corrected Ha luminosity matches that from 
attenuation-corrected UV luminosity (e.g., Buat et al. 2002; Treyer et al. 2007). Using this 
method, Buat et al. (2002) explored the possibility of using the observed Ha/UV ratio to 
estimate the UV attenuation. Treyer et al. (2007) employed this method to estimate the UV 
attenuation and then built an Afuv ~~ FUV-NUV relation. The alternative approach is to 
obtain a scale factor between the attenuations in the Ha line and in the UV by assuming a 
dust attenuation law (e.g., Buat et al. 2002). Both methods are strongly model-dependent, 
however, which adds to systematic uncertainties in the results. The UV attenuation derived 
using the first method, i.e., by matching the UV and Ha derived SFRs, depends not only 
on the Balmer-based attenuation estimates in the Ha line but also on the conversions of Ha 
and UV luminosities to SFRs. On the other hand, the second method rests heavily on the 
form of the assumed dust attenuation law. 

In this paper we take advantage of the large sets of multi-wavelength observations of 
nearby galaxies, including UV, IR and optical spectra derived Balmer decrement ratios, to 
calibrate IR/UV and FUV-NUV as UV attenuation estimators empirically in a self-consistent 
way. The calibrations derived this way are purely empirical, and depend only weakly on the 
Balmer-derived A Hq . 

This paper is the second of a series, in which we aim to use multi-wavelength mea- 
surements of nearby galaxies to derive attenuation corrections. In Kennicutt et al. (2009; 
hereafter Paper I), we used the Spitzer Infrared Nearby Galaxies Survey (SINGS; Kenni- 
cutt et al. 2003) and an integrated spectral survey sample from Moustakas & Kennicutt 
(2006; hereafter MK06) to calibrate the combinations of the optical emission lines (Ha and 
[OII]A3727) and the IR or radio continuum emission as attenuation-corrected SFR indicators. 
We found that the dust attenuation corrected Ha luminosities, as derived from the Balmer 
decrement ratios, can be closely traced by the combinations of the observed Ha emission-line 
luminosities with 24/im, total IR (TIR; the bolometric luminosity over the wavelength range 
3-1100/im), and even 8/zm IR luminosities, with rms scatters ~ 0.1 dex. 

In this paper we derive new calibrations for UV attenuation corrections as derived from 
the IR/UV ratio and the FUV-NUV color as described above, and test the consistency of the 
SFRs computed from IR/UV or FUV-NUV corrected UV luminosities to those from Balmer- 
corrected Ha luminosities. Compared to FUV, NUV can be contributed by more evolved 
stars (~ 1 Gyr). Therefore we mainly focus on FUV in this paper. But we also explore 
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the reliability of using IR/NUV as a dust attenuation estimator as NUV is sometimes used 
as a SFR indicator. The remainder of this paper is organized as follows. In Section 2 we 
describe the sample, and describe how the UV flux densities were obtained. The calibrations 
of the relationships between the attenuation in the FUV and the TIR to FUV ratio, the UV 
color are shown in Section 3. In Section 4 we apply the calibrations derived in Section 3 
to two samples of nearby galaxies to test the reliability of using these attenuation-corrected 
FUV luminosities as SFR indicators. Second order effects of these attenuation-corrected 
FUV luminosities relative to the Balmer-corrected Ha luminosities are also investigated 
in Section 4. At last, we explore the possibility of using the combinations of FUV with 
monochromatic MIR luminosities 25/zm (or 24/xm) and radio continuum to estimate dust 
attenuations in Section 4. In Section 5, the reliability of using NUV luminosity as a SFR 
indicator is discussed. In Section 6 the limitations and the range of applicability of the use of 
the calibrations are discussed. Comparisons of our calibrations to those in the literature, to 
the attenuation curve of Calzetti et al. 2000 and to results based on SFR matching method 
are made in Section 6 as well. Finally we summarize our results in Section 7. For distance 
estimates, we assume H = 70km s -1 Mpc -1 with local flow corrections. 



2. SAMPLE AND DATA 

Our galaxies are drawn from the integrated spectrophotometric survey by MK06 and the 
SINGS survey. The survey information and the selection criteria were described in Paper I. 
We briefly list them here for completeness. The MK06 survey observed 417 galaxies, which 
were selected to cover the full range of optical spectral characteristics found in the local 
galaxies, using a long-slit drift-scanning technique. The reduced spectrum was integrated 
over a rectangular aperture that includes most or all of the emission of the galaxy. The SINGS 
survey was a mult i- wavelength survey from the FUV to the FIR and it was composed of 75 
nearby galaxies that cover wide ranges in morphological type, luminosity, star formation rate 
and dust opacity. In this paper series we aim to use attenuation-corrected Ha luminosities, as 
derived from Balmer decrement ratios, to calibrate combinations of optical emission lines or 
UV continuum and the IR or radio continuum as dust attenuation corrected SFR indicators. 
Therefore several criteria were employed to ensure the quality of the data (1) galaxies without 
detectable Ha emission were excluded, (2) in order to obtain reliable measurements of the 
Balmer decrement ratios, galaxies with low signal to noise ratio spectra at H/3 (S/N < 
15) were removed, (3) galaxies classified as active galactic nuclei (AGNs) or AGN-starburst 
composites (c.f. Moustakas et al. 2006; Kewley et al. 2001; Kauffmann et al. 2003) were 
excluded, and (4) galaxies with lower than 3<x detections at 25/im were removed. Given 
that we will calibrate the attenuation estimates in the UV from the TIR/UV ratio and the 
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FUV-NUV color, FUV and NUV photometric data are needed. Some of the galaxies used 
in Paper I have not been observed by the Galaxy Evolution Explorer (GALEX; Martin et 
al. 2005) so they are not included in the following analyses. 

Similar to Paper I, IR data and the attenuation corrected Ha luminosities from Ha/H/3 
ratios are required. The compilation of Ha fluxes, [Nil] /Ha, Ha/H/3 and IR data has 
been described in detail in that paper. As stated in Paper I, most of the galaxies in the 
MK06 sample have not been observed with Spitzer. To maximize the uniformity of the 
data, only the Infrared Astronomical Satellite (IRAS) data were used to derive the TIR 
luminosities 1 , which was estimated using a weighted sum of IRAS 25yum , 60/im and 100/im 
emission following Dale & Helou (2002; equation (5) of their paper). We will discuss possible 
systematics introduced by using IRAS-based TIR luminosities in Section 6.6. We also provide 
calibrations for the combinations of the monochromatic MIR (IRAS 25/iin or Spitzer 24/zm) 
luminosities and the observed UV luminosities in Section 4.3. This is useful when FIR data 
are not available. As shown in Paper I, the IRAS 25/im luminosity (L(A) = zaL>) and the 
Spitzer/Multiband Imaging Photometer for Spitzer (MIPS) 24/im luminosity can be used 
interchangeably. Given that Spitzer/MIPS 24/im photometry is more sensitive than IRAS 
we use those measurements in preference to IRAS for the SINGS galaxies in Section 4.3. 
Spitzer/MIPS 24/im measurements were also performed for the central 20" x 20" regions of 
the SINGS galaxies in Paper I, and they are also used in this paper (Section 4.3). 

For galaxies selected from the MK06 sample, we identified their counterparts in the 
UV from GALEX GR4. GALEX performs imaging observations at FUV and NUV bands, 
centered at 1528A and 227lA respectively, with a circular field of view of 1?2 diameter 
(Martin et al. 2005). GR4 includes images taken from the All-sky Imaging Survey (AIS), 
Medium Imaging Survey (MIS), Nearby Galaxies Survey (NGS), Deep Imaging Survey (DIS) 
and Guest Investigators Survey (GI). Since artifacts appear much more frequently at the edge 
of the field of view, only galaxies within a radius of 0?55 from the field center were selected 
(Morrissey et al. 2007). Because of the overlap in coverage of these surveys, sometimes 
there were more than one image available. In such cases, we used the one with the longest 
exposure time in the FUV. After this process we were left with 98 galaxies (54 in AIS, 4 in 
MIS, 22 in NGS, 2 in DIS and 16 in GI) out of the 113 star-forming galaxies selected from 
the MK06 sample. But NGC 1003, which is in AIS, was excluded later because of artifact 
contamination in its FUV image. So the final sample consists of 97 objects, with 53 in the 
AIS survey. 



1 As shown in Paper I, the TIR luminosities derived using IRAS measurements are systematically lower 
than those based on Spitzer MIPS measurements by 24% ± 23% on average. 
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Intensity (-int) and sky background (-skybg) images at FUV and NUV bands were re- 
trieved from the Multimission Archive at Space Telescope Science Institute (MAST) website 2 
(Morrissey et al. 2007). Sky background was subtracted before performing photometry. In 
order to obtain the UV fluxes from the same region as Ha, the rectangular apertures used in 
the integrated spectra observations by MK06 were adopted in the UV photometry 3 . Fore- 
ground stars within the apertures, identified by comparing the NUV and the Digitized Sky 
Survey (DSS) optical images, were masked during the photometry. The FUV and NUV 
magnitude zero points listed in Morrissey et al. (2007) were used to convert intensities into 
magnitudes. The high-resolution relative response images (-rrhr) from MAST were used to 
estimate the uncertainties in the photometry. The FUV and NUV flux densities were cor- 
rected for Galactic extinction using the Schlegel et al. (1998) dust map and the Galactic ex- 
tinction curve derived by Cardelli et al. (1989) for a total-to-selective extinction of Ry = 3.1 
(c.f., Gil de Paz et al. 2007). Specifically, A FUV = 7.9E(B - V) and A NUV = 8.0E(B - V). 
The Galactic extinction corrected FUV and NUV flux densities of the MK06 sample are 
listed in Table 1. The absolute calibration uncertainties in FUV (0.05 mag) and NUV (0.03 
mag) are added in quadrature. 

Since more than half of the MK06 sample galaxies (53 out of 97) were extracted from 
the shallow AIS survey, with exposure times of 100s, 15 times shorter than the MIS, we 
tested whether some of the extended UV emission was missed in AIS images compared to 
observations with deeper exposures. There are 43 objects in common between the AIS and 
the other deeper surveys. We compared their UV magnitudes measured from the AIS images 
and those from observations with longer exposure time in Figure 1. It is clear that the UV 
measurements based on the AIS images are in excellent agreement, within ~20%, with those 
from deeper exposures, which means that the UV fluxes measured from the AIS images are 
not systematically under-estimated within the adopted rectangular apertures. Consequently 
they can be used along with the measures from deeper exposures without introducing any 
systematic errors. 

UV flux densities for the SINGS sample were obtained from Dale et al. (2007). Those 



2 http: / /galex.stsci.edu/ GR4/ 

3 We note that the IR fluxes we used were measured by IRAS and hence they are not aperture-matched. 
Because of the poor resolution of IRAS imaging, we cannot estimate the fraction of IR emission outside of the 
rectangular apertures from the IRAS imaging data. Instead we used Spitzer/MIPS 24fim images for galaxies 
in both MK06 and SINGS samples to get a rough estimate of this fraction. By comparing the rectangular 
photometry to the global photometry which was adopted from Dale et al. (2007), we found that at most 
10% of 24/um emission is from outside of the rectangular apertures. This small fraction will not affect our 
results. 
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authors calculated the FUV and NUV flux densities from the surface brightness profiles de- 
rived for the GALEX Atlas of Nearby Galaxies (Gil de Paz et al. 2007). Since Dale et al. 
(2007) used Li & Draine (2001) reddening curve to correct for the Galactic extinction, we 
redid the Galactic extinction using the same recipes as those used for the MK06 sample for 
consistency 4 . Applying the requirement of having both IRAS data and UV flux measure- 
ments left us with 36 objects for the SINGS sample 5 . To avoid small number statistics and 
the less reliable measurements of Ha/H/3 ratios for some of the SINGS objects (see Paper 
I), we do not use these objects to calibrate the relations in the following analyses. However 
we include them in the comparison to see how they behave. We also measure the FUV and 
NUV flux densities for the central 20" x 20" regions of the SINGS galaxies. It should be 
noted that the large absolute calibration uncertainties (0.15 mag) in Gil de Paz et al. are 
included in the quoted errors in FUV and NUV for SINGS sample. 



3. EMPIRICALLY CALIBRATED A FUV - TIR/FUV, FUV-NUV COLOR 

RELATIONS 

As addressed in Section 1, prescriptions of using TIR/FUV as a dust attenuation in- 
dicator are usually derived under various model assumptions. In this section, we employ 
a new method to derive such prescriptions empirically and avoid most of the systematic 
model dependencies. Our method is based on the statistical correlation between the FUV 
attenuation estimated from the TIR/FUV ratio and that as derived from the FUV-NUV 
color. By deriving the prescriptions based on a sample of nearby star-forming galaxies, we 
virtually assume an average stellar population and dust attenuation curve of a certain type 
of galaxies. Therefore our calibrations are probably not applicable for other types of galaxies 
(see Section 6). 

We adopt the convention of infrared excess IRX from Meurer et al. (1995) 

IRX = log[L(TIR)/L(FUV) obs ], (1) 
where the observed FUV luminosity takes the common definition of monochromatic lumi- 



4 Due to a mistake in the code, the UV flux densities for the SINGS sample published in Dale et al. (2007) 
were virtually foreground corrected according to what Gil de Paz et al. (2007) used for a reddening curve, 
same as those used for the MK06 sample. So they were used directly to compute luminosities and colors. 
No corrections were made. 

5 Since MIPS 24/im measurements are used in preference to IRAS 25/im for the combination of FUV and 
25yum (or 24/xm) for the SINGS galaxies in Section 4.3, all star-forming galaxies in the SINGS sample (see 
Paper I) with GALEX observations are included. 47 objects meet this requirement. 
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nosity L(FUV) obs = i/L„(1528A). 

According to an energy balance argument, the dust attenuation in the FUV, A FUV , can 
be estimated from the IRX using the following form 

A FUV = 2.5 log(l + a FUV • 10 IRX ), (2) 

where a F uv is a scale parameter, whose meaning is explained below. This approach is similar 
to that presented in Paper I for the combination of the observed Ha luminosity and the TIR 
luminosity. As stated in Paper I, such an expression is a simple approximation to a much 
more complicated radiative transfer process. To better understand the underlying physical 
foundations of eq. (2), we adopt the relevant equations from Paper I and adapt them to the 
case of FUV. 

We first define a scale parameter rj, which is the fraction of the bolometric luminosity 
that is emitted in the FUV, i.e., the inverse of the bolometric correction 

L(FUV) corr = ?7FuvLboi, (3) 

where L(FUV) corr represents the intrinsic (extinction-free) FUV luminosity and L bo i denotes 
the bolometric luminosity. 

Under the assumption of a foreground dust screen approximation, the observed lumi- 
nosity can be expressed as 

L(FUV) obs = L(FUV) corr e-^ uv , (4) 

where tfuv is an effective optical depth in the FUV, as defined by equation (4), and the 
corresponding attenuation in magnitudes Apuv = 1.086 tfuv- 

The dust attenuation corrected FUV luminosity is a sum of the observed luminosity and 
the dust-attenuated luminosity 

L(FUV) corr = L(FUV) obs + L(FUV) corr (1 - e"**™). (5) 



Considering the energy conservation, all the starlight reprocessed by dust is re-emitted 
in the IR. So the TIR luminosity can be expressed as 

L(TIR) = L bol (l-e- f ), (6) 

where f is the effective opacity for all of the dust-heating starlight, and is defined as 

/ L x e- T *d\ 
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This effective mean opacity depends on the star formation history of the galaxy, the dust 
extinction curve and the relative distribution of stars and dust. 

We substitute eqs. (3) and (6) into eq. (5) to eliminate L bol and express L(FUV) corr us- 
ing the observable luminosities L(FUV) obs and L(TIR) and the opacities. After this exercise, 
the dust attenuation estimate in the FUV becomes 

A FUV = 2.5 log[l + 7? FU v • 10 (1 _ e ^ ]• ( 8 ) 

By comparing eq. (2) with eq. (8), it can be seen that the coefficient ofuv in eq. (2) is 
the approximation of the product of the inverse of the bolometric correction ?7fuv and the 
factor ( Yi_ e T y • 

When the FUV-NUV color is used to estimate the FUV attenuation, a linear relation 
between A F uv and the FUV-NUV color is usually taken 

Afuv = ^fuv ■ (FUV - NUV) obs + (9) 

where the slope sfuv and the intercept 2fuv are usually determined by a linear regression 
to the observations. The physical foundations of this relation can be understood as follows. 

The color excess in FUV-NUV, E(FUV-NUV) is defined as 

E(FUV - NUV) = (FUV - NUV) obs - (FUV - NUV) int = A FUV - A NUV , (10) 

where (FUV — NUV) obs and (FUV — NUV) int denote the observed and the attenuation- 
corrected FUV-NUV colors. E(FUV-NUV) relates to the color excess E(B-V) by the relation 

E(B-V) = - — -E(FUV-NUV), (11) 

(kFUV — kNuv) 

where kA = Aa/E(B — V) represents an attenuation curve. So eq. (9) has the modified form 
Afuv = T , ^ r[(FUV - NUV) obs - (FUV - NUV) int ]. (12) 

(^FUV — -KNUV J 

By comparing this modified equation to the original eq. (9), we can see that sfuv = 
(k FU t F -k NUV ) ' w hich is the slope of the UV part of the attenuation curve, and ?fuv/sfuv 
is equal to the dust-free UV color, (FUV - NUV) int . 

Since the estimated Apuv from IRX and FUV-NUV should agree, the predicted relation 
between IRX and FUV-NUV color can be obtained by equating eq. (2) and eq. (12). Given 
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that kpuv and Icnuv cannot be derived separately from the fitting, sfuv is still used in the 
1RX - FUV-NUV relation 

1RX = l og [10°- 4s ™v-[(FUV-NUV) obs -(FUV-NUV) int ] _ ^ _ log(aFUV ). ( 13 ) 

Before we proceed to derive the parameters ofuv, (FUV — NUV)i nt and sfuv in eq. (13) 
from the observations, we discuss the implicit assumptions in our method. The comparison 
between eq. (2) and eq. (8) shows that the coefficient ofuv is the approximation of the 
product of the inverse of the bolometric correction r^puv and the factor ^l-e^)^ • By 
deriving single values for apuv; (FUV — NUV) int and spuv, we virtually study the average 
properties of this sample, without taking into account variations in star formation histories 
and dust attenuation curves from galaxy to galaxy, which contribute to the scatters in the 
observed correlations. 

In order to minimize the number of free parameters in the fitting, we first estimate the 
intrinsic (i.e., dust-free) FUV-NUV color empirically, which is achieved by the correlation 
between the observed FUV-NUV color and the attenuation in the Ha line. This method 
is similar to Calzetti et al. (1994) in which the unattenuated power-law slope of the UV 
continuum /3 was estimated using the correlation between j3 and the Balmer decrement 
ratio. This correlation between the dust reddening in the stellar continuum and in the gas 
content (see also Calzetti 1997 and Figure 2 below) implies that the FUV attenuation in 
the stellar continuum is proportional to the attenuation in the Ha line. The proportional 
factor is dependent on the extinction curve and the relative distribution of stars, gas and 
dust, which may vary from galaxy to galaxy. Unfortunately, we cannot study this property 
on a galaxy-to-galaxy basis based on current data, so only an average relation is obtained 
(see Section 6.3). 

The FUV-NUV color versus Ah q relation is shown in Figure 2. Throughout this paper 
the MK06 and SINGS objects are represented by solid and open circles respectively. From 
this figure, we can see that the MK06 and the SINGS samples tell us different stories. The 
MK06 sample shows that Ah« is correlated with FUV-NUV color, whereas the SINGS sample 
shows almost no correlation between these two quantities after the two reddest objects are 
excluded. This is not unexpected given the large calibration uncertainties of UV photometry 
and the way we compiled the Balmer decrement ratios for the SINGS galaxies, as described 
in Section 2 and Paper I. Therefore only the MK06 sample is used to derive the dust-free 
FUV-NUV color, as indicated in Section 2. Since the scatter is probably intrinsic, which 
means that the measurement errors are not the main contributors to the scatter about the 
regression line, an ordinary least square bisector fit method is adopted (Isobe et al. 1990; 
Feigelson & Babu 1992). The fitting to the MK06 sample that is plotted as a solid line in 
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Figure 2 gives 



(FUV - NUV) = (0.539 ± 0.029)A Ha + (0.022 ± 0.024) 



(14) 



This indicates an intrinsic FUV-NUV color for zero attenuation of 0.022 ± 0.024 mag. The 
errors quoted here in the slope and intercept were yielded by the fitting routine. The dust- 
free FUV-NUV color derived here is in excellent agreement with the value 0.025 ± 0.049 
mag derived by Gil de Paz et al. (2007). Those authors used the prescriptions given by 
Buat et al. (2005) to correct for the dust attenuation in the FUV and NUV using the 
TIR/FUV and TIR/NUV ratios, and found that the intrinsic FUV-NUV color for their 
sample is 0.025 ±0.049 mag on average. To understand the average stellar population of our 
sample galaxies, we compare the intrinsic FUV-NUV color we derived here with the stellar 
population synthesis predictions. By comparing to models constructed using GALAXEV by 
Bruzual & Chariot (2003), we find that our galaxies are ~ 0.02, 0.03 and 0.04 mag redder in 
FUV-NUV color than a 10 Gyr old stellar population with an e- folding timescale of 5 Gyr, 10 
Gyr and a constant star formation history, respectively. For a constant star formation history 
STARBURST99 (Leitherer et al. 1999; Vazques & Leitherer 2005) produces consistent result 
(see Section 6.6). 

After fixing the term (FUV — NUV)i nt = 0.022 in eq. (13), the coefficients sfuv and 
cifuv can be obtained by fitting eq. (13) to the data. Figure 3 shows the IRX versus 
FUV-NUV color relation, with the corresponding /3 labelled on the top axis. The relation 
between FUV-NUV and /3 was derived using the definition given by Kong et al. (2004) for 
our adopted central wavelengths for FUV and NUV (see Section 2), which has the form 
j3 = 2.32(FUV — NUV) — 2.00. In order to have a clearer view of the average properties of 
the data, we overplot the IRX median, lower (25%) and upper (75%) quartiles in bins of 
width 0.2 mag in FUV-NUV color as open squares with error bars in Figure 3. The solid line 
represents the chi-square minimization fit to the MK06 sample with the deviation in the IRX 
minimized, which has the form defined by eq. (13) with sfuv and cifuv equal to 3.83 ± 0.48 
and 0.46±0.12 respectively. In order to test the effect of the adopted (FUV — NUV)i nt on the 
fitting results, we varied (FUV — NUV) int by ±la about the best-fitted value, i.e., assumed 
(FUV - NUV) int = 0.046 and (FUV - NUV) int = -0.002. It turned out that the changes in 
the best fitted sfuv and a F uv caused by these new assumptions of the (FUV — NUV) int are 
within ±1<7, and the resulting IRX versus FUV-NUV color relations cannot be distinguished 
from that using (FUV - NUV) int = 0.022 over the ranges of FUV-NUV color and IRX 
covered by our sample galaxies (see Section 6.1). Substituting the best fitted coefficients 
into eqs. (2) and (12), the attenuation estimates based on the IRX and the FUV-NUV color 
become 



A 



FUV — 



2.5 log[l + (0.46 ± 0. 12) • 10 mx ] , 



(15) 
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and 

A FUV = (3.83 ± 0.48) [FUV - NUV) obs - (0.022 ± 0.024)]. (16) 

4. DUST ATTENUATION CORRECTED FUV LUMINOSITY AS A SFR 

INDICATOR 

In order to evaluate the reliability of the dust attenuation corrected FUV luminosities 
derived using the above calibrations as SFR indicators, we compare the IRX and FUV-NUV 
color corrected FUV luminosities (according to eq. (15) and (16)) with the Balmer decrement 
ratio corrected Ha luminosities and investigate possible second-order trends in the residuals 
as functions of various properties of the galaxies in this section. Combinations of FUV 
with monochromatic 25/im (or 24yum) infrared and radio continuum as dust attenuation 
measures are also explored. We note that the attenuation calibration methods we employed 
in Section 3 make only minimal use of the Balmer decrement data. So when we compare the 
dust-corrected FUV luminosities based on our calibrations with the Balmer-corrected Ha 
luminosities, we are not making a circular argument. 

Before we proceed to compare the dust attenuation corrected FUV luminosities with 
the Ha luminosities, it is interesting to see the behavior of the uncorrected FUV luminosities 
as a function of the Balmer decrement corrected Ha luminosities in Figure 4, which can help 
us understand the importance of the attenuation corrections for this sample. The straight 
lines in Figure 4 are model predictions for matched SFRs as described below. Figure 4 shows 
that the uncorrected FUV luminosities can underestimate the SFRs by more than an order 
of magnitude, which varies with the intensity of the star formation activity. Galaxies with 
different SFRs are not equally attenuated, with more active galaxies suffering from higher 
attenuations. This is consistent with studies in the literature (e.g., Wang & Heckman 1996; 
Treyer et al. 2007). In the following part of this section, the reader will see how much will 
be improved with our dust attenuation corrections. 

4.1. Infrared-Corrected FUV Luminosities As SFR Indicators 

We compare the IRX-corrected FUV luminosities with the Balmer decrement ratio cor- 
rected Ha luminosities in the left panel of Figure 5. It can be clearly seen that the IRX- 
corrected FUV luminosities correlate tightly and linearly with the Balmer-corrected Ha 
luminosities with a rms scatter of ±0.09 dex. 

In order to quantify the degree of consistency between the SFRs derived from the 
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attenuation-corrected FUV and Ha luminosities, in the left panel of Figure 5 we overplot 
the predicted relations from different SFR prescriptions. The red short-long dashed line 
represents the prediction from the widely used K98's SFR prescriptions. The calibrations 
in K98 assume a Salpeter initial mass function (IMF; Salpeter 1955) with mass limits 0.1- 
100M Q , and 1990's generation stellar evolution models (see Kennicutt et al. 1994; Madau et 
al. 1998). To compare with the K98's prediction, we first construct a model under the same 
assumptions as K98 did, i.e., Salpeter IMF, constant star formation history lasting for 100 
Myr and solar metallicity (the magenta short dashed line), using Version 5.1 of the stellar 
population synthesis model STARBURST99 (Leitherer et al. 1999; Vazquez & Leitherer 
2005), which uses the state-of-the-art stellar evolutionary models. 

Star-forming galaxies in our sample probably have experienced continuous star forma- 
tion for at least 1 Gyr. In order to examine age effects, we build models with a constant star 
formation history lasting for 1 Gyr (represented by long dashed lines). It is worth noting 
that for a constant star formation history the intrinsic FUV to Ha luminosity ratio stays 
constant after 1 Gyr. Therefore, the model predictions we derive by assuming a constant 
SFR for 1 Gyr are the same as those under the assumption of a constant SFR for longer 
time. Although we assume constant star formation histories here, models with exponentially 
declining star formation histories do not give significantly different results, as long as the 
SFR does not drop too rapidly. For example, for an e-folding timescale of 5 Gyr, the FUV 
to Ha ratio differs from the constant star formation history by < 4% at any given age, 
according to our experiments based on stellar population synthesis models constructed using 
GALAXEV by Bruzual & Chariot (2003). 

The IMF also affects the FUV/Ha ratio. So apart from models with Salpeter IMF, 
models with a more realistic Kroupa IMF (Kroupa & Weidner 2003) with mass limits of 
0.1-IOOMq are also explored (represented by lines in blue). The FUV and Ha luminosities 6 
and their ratios are listed in Table 2. 

In the left panel of Figure 5 the different models (represented by the colored lines) are 
difficult to separate, because of the large dynamic range in luminosities. For that reason 
we show in the right panel of Figure 5 histograms of the attenuation-corrected FUV to Ha 
luminosity ratio, using a logarithmic scale. The MK06 and SINGS samples are represented 
by the solid and dotted histograms respectively. The straight lines in color are the same 
model predictions as those shown in the left panel of this figure. The median and mean 
values of the IRX-corrected FUV to Ha luminosity ratios for the MK06 sample are 123.88 



6 In fact we do not list the luminosities themselves. Instead we give the scale parameters - Cfuv and 
Cua, which connect luminosities and their corresponding SFRs by the relation SFR(A) = C\ ■ L(A). 
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and 127.94, which are in excellent agreement with the model predictions for Kroupa IMF at 
100 Myr and Salpeter IMF at 100 Myr respectively (see Table 2). However, the FUV/Ha 
ratio only increases by 8% and 9% respectively for the same IMFs but older age (> 1 Gyr). 

The correlation between the IRX-corrected FUV luminosities and the Balmer-corrected 
Ha luminosities is so tight and linear that it is useful to test for any second order trends 
in the residuals as functions of various galactic properties. In Figure 6, we plot the ratios 
of the IRX-corrected FUV luminosities to the Balmer-corrected Ha luminosities against the 
Balmer-corrected Ha luminosities (i.e., the dust corrected SFRs; top- left panel), the Balmer- 
corrected Ha luminosity densities (top-right panel), the attenuations in the Ha line (bottom- 
left) and the axial ratios b/a (i.e., inclinations; bottom-right). In Figure 7, the optical 
spectral features - Ha equivalent width EW(Ha) (top-left panel), 4000A break D n (4000A) 
(top-right panel), gas-phase oxygen abundance 12+log(0/H) (bottom-left panel) and the 
dust temperature indicator - far-infrared 60/im to 100/im flux ratio / I/ (60/im)// I/ (100/im) 
(bottom-right) are explored. Given that Ha line and continuum emission are SFR and stellar 
mass tracers respectively, the integrated Ha equivalent width is a measure of specific SFR 
(Kennicutt et al. 1994). The D n (4000A) is often used reddening-insensitive rough star 
formation history indicator of the stellar populations (e.g., Bruzual 1983; Balogh et al. 1999; 
Kauffmann et al. 2003) and it correlates with the specific SFR (Brinchmann et al. 2004). In 
all panels the dashed line indicates the predicted FUV to Ha luminosity ratio for a Kroupa 
IMF at age 100 Myr, as shown in Figure 5. 

As can be seen from Figure 6 and 7, the ratios of the IRX-corrected FUV to the 
attenuation-corrected Ha luminosities show no or at most marginal trends with all the pa- 
rameters investigated. The rms scatters are about ±0.09 dex, which is even smaller than 
the average uncertainty in log [L(FUV)iRx,corr/L(Ha) corr ] ratio, ±0.13 dex, with the error 
introduced by using eq. (15) to correct for FUV attenuation included. This demonstrates 
the reliability of using the combination of the FUV and TIR luminosities as a SFR indicator. 
The absence of correlations implies that over the ranges of the galactic properties probed by 
our galaxies, either no trend exists in the dust-free FUV/Ha ratios and the IRX method with 
the parameters tested above, or any trend in the dust-free FUV/Ha ratios is compensated 
by that in the IRX method. Whatever it is, the SFR estimated from the IRX-corrected FUV 
luminosity agrees well with that derived from the Balmer-corrected Ha luminosity, and the 
ratio of these dust-corrected SFRs is not a function of the galactic properties. Therefore the 
combined FUV and TIR luminosities can serve as a very good attenuation-corrected SFR 
indicator over the ranges of physical properties covered by our sample objects. 

From Figure 5 to Figure 7, we note that the SINGS galaxies deviate from the MK06 
galaxies systematically. The systematic deviations can be caused either by overestimates of 
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the TIR luminosities or by underestimates of the attenuation in Ha lines. As we showed in 
Paper I (Figure 2 in that paper; see also Dale et al. 2009), for cold galaxies (/ 1/ (60/xm)// 1/ (100/im) < 
0.6), the TIR luminosities measured by IRAS bands are systematically underestimated with 
respect to the more sensitive and accurate Spitzer/MIPS bands. Therefore using TIR lu- 
minosities based on Spitzer/MIPS measurements for the SINGS galaxies would make the 
SINGS galaxies deviated more. Furthermore, we did not see a correlation between the devi- 
ations and the / I/ (60/im)// J/ (100/im) colors in the bottom-right panel of Figure 7. Another 
piece of evidence is that the same deviations show up in the comparison of the combined 
FUV+25/im luminosities with the Balmer-corrected Ha luminosities (Section 4.3), and the 
monochromatic 25/xm luminosity measured by IRAS does not suffer from the same problem 
as the TIR luminosity 7 . Therefore, the systematic deviations of the SINGS galaxies cannot 
result from overestimates of the TIR luminosities by IRAS. Instead, they may be due to 
underestimates of the attenuation in Ha lines, caused by the noisier H/3 detections. 



4.2. FUV-NUV Corrected FUV Luminosities As SFR Indicators 

Similar to the analyses we did for the IRX-corrected FUV luminosities in Section 4.1, 
we compare the FUV-NUV color corrected FUV luminosities with the Balmer decrement 
ratio corrected Ha luminosities in this subsection. Figure 8 shows the FUV-NUV corrected 
FUV luminosities as a function of the Balmer-corrected Ha luminosities (left panel) and 
the histogram of the FUV-NUV corrected FUV to the Balmer-corrected Ha luminosity ratio 
(right panel). Compared to the IRX-corrected FUV luminosities versus the Balmer-corrected 
Ha luminosities relation, the correlation between the FUV-NUV corrected FUV luminosities 
and the Balmer-corrected Ha luminosities is non-linear and has a much larger rms scatter 
of ±0.23 dex. 

Figure 8 shows that the relation between the FUV-NUV color corrected FUV luminosi- 
ties and the Balmer-corrected Ha luminosities is not completely linear. The slope of the 
relation is less than unity, which implies that for galaxies forming stars vigorously (i.e., with 
high SFRs), the FUV-NUV corrected FUV luminosities will under-estimate the SFRs. 

To better understand the applicability and limitation of the FUV-NUV corrected FUV 
luminosities as SFR measures, apart from the Balmer-corrected Ha luminosities, we also ex- 
amined the ratios of the UV color corrected FUV luminosities and the Balmer-corrected Ha 
luminosities as functions of several other galactic properties, as we did for the IRX-corrected 



7 In paper I, we found that the IRAS 25^m luminosities are in well agreement with the Spitzer/MIPS 
24/xm luminosities. 
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FUV luminosities. We found that the residuals of the UV color corrected FUV luminosi- 
ties relative to the Balmer-corrected Ha luminosities 8 show significant negative correlations 
with several of the parameters, including the Balmer-corrected Ha luminosity densities, the 
attenuations in the Ha lines, the FIR color and the EW(Ha), as shown in Figure 9 and 10. 
The negative trends with the FIR color and the Ha equivalent width are similar to what 
Kong et al. (2004) found. Those authors found correlations between the deviation from 
the average IRX-/3 relation defined by starburst galaxies and the above two parameters, and 
they explained these correlations as direct observational evidence for the dependence of the 
scatter of the IRX-/3 relation on star formation history. Galaxies that are forming stars more 
actively in the current epoch than in the past are intrinsically bluer than those with lower 
current to past SFR ratios. So the attenuation in those galaxies with warmer FIR color, 
larger EW(Ha) and hence bluer intrinsic FUV-NUV color, would tend to be under-estimated 
by the prescription defined by the average population of the galaxy sample (see eq. (16)). 
Since stars older than ~ 1 Gyr do not contribute to UV emission, the above argument is only 
valid for galaxies with either continuous star formation history or continuous star formation 
history with superposed bursts younger than ~ 1 Gyr. For a continuous star formation 
history with superposed bursts taking place at more than ~ 1 Gyr ago, the FUV-NUV color 
behaves as if no bursts occurred at all in the history. On the other hand, the correlations 
of the residuals of the FUV-NUV corrected FUV luminosities with the Ha luminosities, the 
Ha luminosity densities and the Ha attenuations are hard to be interpreted in terms of 
the dust-free FUV-NUV color because there is no evidence for a correlation between the 
absolute star formation activity and the intrinsic FUV-NUV color. The possible reason for 
these correlations is the change in the effective attenuation curves (see eq. (12)) with the 
star formation intensities given that the intensity of the star formation activity may change 
the dust properties and the stars/dust geometry, and hence the effective attenuation curves. 
These speculations agree with the argument by Burgarella et al. (2005) and Boquien et al. 
(2009), who suggested that both star formation history and dust attenuation curve effects 
should be taken into account when the FUV-NUV color is used as a UV attenuation esti- 
mator. In other words, a simple formalism of using FUV-NUV color as a FUV attenuation 
measure (e.g., eq. (9)) is far from enough. Additional constraints on star formation histories 
and dust attenuation curves are needed. 



8 In order to minimize any systematic effects introduced by using the Balmer-corrected Ha luminosities 
as references, we also examined the trends of the residuals of the UV color corrected FUV luminosities 
relative to the IRX-corrected FUV luminosities as functions of the above parameters. We found that these 
residuals behave similarly to the residuals of the UV color corrected FUV luminosities with respect to the 
Balmer-corrected Ha luminosities. 
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4.3. Other Composite SFR Indicators 

Since FIR photometry and hence the TIR luminosity is sometimes not available, it is 
worthwhile calibrating the combinations of the FUV and the monochromatic IR waveband 
luminosities (e.g., 25/im luminosities 9 ) and radio emission at 1.4GHz as attenuation-corrected 
SFR measures. 

By matching the combination of the observed FUV and the 25/im infrared luminosi- 
ties with the IRX-corrected FUV luminosities, we obtain a coefficient of 3.89 ± 0.15 for 
FUV+25/zm , which indicates an average L(25/im ) to L(TIR) ratio of ~ 0.12. This is well 
within the observed range for star- forming galaxies (c.f., Paper I; Calzetti et al. 2010). The 
left panel of Figure 11 shows a tight and linear correlation between the Balmer-corrected Ha 
luminosities and the combined FUV and 25/zm luminosities. Since the relation between the 
Balmer-corrected Ha luminosity and the IRX-corrected FUV luminosity, which is the refer- 
ence of the combined FUV and 25/zm luminosity, is consistent with the predicted relation 
for a constant star formation history lasting for 100 Myr and Kroupa IMF (see Section 4.1), 
we plot this model prediction here for reference (the dashed line). As expected, the data 
match the model well. The dispersion is ±0.13 dex. 

In this panel, besides SINGS (open circles) and MK06 (filled circles) samples, we also 
include the central 20" x 20" regions of the SINGS galaxies, as denoted by open squares. 
Although the calibration was derived based on the integrated measurements of MK06 sample, 
it also compares well to the centers of the SINGS galaxies. This perhaps is not surprising 
because of the large coverage of the centers (typically several square kilo-parsecs) . 

To our knowledge, there exists only one such calibration in the literature. Zhu et al. 
(2008) calibrated a FUV and 24/zm combination by comparing to a FUV-NUV corrected FUV 
luminosity. They used the prescription published by Treyer et al. (2007) to derive the FUV- 
NUV corrected FUV luminosity and obtained L(FUV) corr = L(FUV) obs + 6.31L(24//m). The 
coefficient of 6.31 is much higher than our value. In order to understand the discrepancy 
between Zhu et al. and our calibration, we compare our eq. (16) with the relation given by 
Treyer et al. We find that the relation provided by Treyer et al. gives systematically lower 
attenuation estimates than ours, which implies that a coefficient lower than 3.89 is required 
to match the FUV-NUV corrected FUV luminosity derived using Treyer et al.'s relation. On 
the contrary, Zhu et al. obtained a much higher value. We cannot figure out the reason that 
leads to a high coefficient obtained by Zhu et al. (2008). 



9 Since the number of galaxies with available data is small, no attempt was made to calibrate the 
combination of UV and 8[im luminosity as an attenuation or a SFR indicator. 



- 19 - 



The combination of the observed FUV and 1.4GHz radio luminosities yields a similarly 
tight relation ±0.14 dex as shown in the right panel of Figure 11, but its correlation with 
the Balmer-corrected Ha luminosities shows clear nonlinearity, with a slope of 1.09 ± 0.02 
from a bisector fit. For the average Ha luminosity of our sample, L(Ha) corr = 10 4L5 L Q , the 
radio-corrected FUV luminosity can deviate from a linear relation by 0.02 dex. We conclude 
that when FIR data are absent, the monochromatic 25/im luminosity and the 1.4GHz radio 
luminosity can be used as a slightly less accurate surrogate of the TIR luminosity to correct 
for the dust attenuation in FUV. 



5. CAN NUV LUMINOSITIES BE USED AS A SFR INDICATOR? 

In the previous sections, we have calibrated the FUV luminosity as a SFR indicator. A 
relevant question is whether NUV can also serve as a reliable SFR tracer. In theory, FUV 
traces emission from stars younger than ~ 300 Myr (e.g., Meurer et al. 2009), whereas 
NUV can be contributed by stars as old as ~ 1 Gyr. FUV and NUV emissions are both 
dominated by stars younger than ~ 100 Myr though. Our experiment with STARBURST99 
shows that for a galaxy with constant star formation history lasting for more than 1 Gyr, 
its FUV emission increases 8% from 100 Myr to 300 Myr old and only 1% from 300 Myr 
to 1 Gyr old. The equivalent numbers for NUV are 12% and 5%, respectively. We note 
that models built from GALAXEV by Bruzual & Chariot (2003) produce consistent results. 
In other words, the NUV emission of such a galaxy contributed by stars older than 100 
Myr is about 17% (compared to 9% for FUV emission), which makes it a less reliable SFR 
tracer in galaxies with a significant fraction of stars being 100 Myr - 1 Gyr old. In the 
past, however, when FUV observations were rarely available NUV measurements were the 
only UV SFR indicators available. Even today, NUV is sometimes used to probe the star 
formation activity (e.g., Iglesias-Paramo et al. 2006); Recipes for using the IR/NUV ratio 
to estimate the attenuation in the NUV have been presented recently (e.g., Buat et al. 2005; 
Burgarella et al. 2005; Cortese et al. 2008). Thus it is useful to assess the reliability of using 
NUV luminosity as a SFR indicator, using our new approach and datasets. 

Given that NUV emission can be from more evolved stars than FUV, we did not fit the 
corresponding coefficient snuv i n eq. (13). Instead we estimated snuv from sfuv according 
to snuv — sfuv — 1 (i-e.,«Nuv — 2.83) 10 . Then cznuv was estimated from fitting eq. (13) 
with NUV substituted for FUV data and snuv fixed. The resulting fitted value of onuv 



10 In fact when we applied the same method as in Section 3 to obtain snuv and (znuv, we got a smaller 
than sfuv — 1 value but this value is still consistent with sfuv — 1 within 1 a error. 
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was 0.27±0.02. Following the analyses for FUV, we compare the TIR/NUV corrected NUV 
luminosities 11 with the Balmer-corrected Ha luminosities in Figure 12. To our pleasant 
surprise, the combined NUV and TIR luminosities correlate with the Balmer-corrected Ha 
luminosities nearly as well as FUV, and the scatter in NUV+TIR vs. L(Ha) corr is almost 
the same as that in FUV+TIR vs. L(Ha) corr (0.10 vs. 0.09 dex respectively). The average 
ratio of the TIR/NUV-corrected NUV luminosity to the Balmer-corrected Ha luminosity 
is consistent with that expected for a constant star formation history over > 1 Gyr for a 
Kroupa IMF. 

We also examined the trends of the residuals of the combined NUV and TIR luminosities 
relative to the Balmer decrement ratio corrected Ha luminosities as functions of various 
parameters. We find that all the behaviors shown in NUV are similar to those in FUV. 
When 25/iin infrared and radio continuum 1.4GHz luminosities are used as substitutes for 
TIR luminosities, we obtain coefficients a NUV of 2.26±0.09 and 41.75±2.97x 10 20 respectively. 
The scatters shown in NUV+25/im and NUV+1.4GHz versus L(Ha) corr are identical to those 
shown in the same comparisons for FUV. 

From the above analyses we can see that the NUV luminosity can serve as a rather 
good SFR indicator for our sample galaxies, and star formation activities have been going 
on smoothly in these galaxies over the last 1 Gyr. However, the use of NUV-based SFRs in 
galaxies with non-constant star formation histories during the last 1 Gyr may be problematic. 
Our experiment with the stellar population synthesis modelling code GALAXEV (Bruzual & 
Chariot 2003) shows that for a galaxy with a constant star formation history lasting for 100 
Myr, the NUV emission at 300 Myr is 4.3% of its value at 100 Myr (compared to 2.6% for 
FUV). For a galaxy with high SFR (~ 23 M yr" 1 ) lasting for 100 Myr, the NUV emission 
at 300 Myr cannot be distinguished from a Milky Way like star-forming galaxy (with SFR 
~ 1 Meyr- 1 ). 



6. DISCUSSION 

In previous sections we have derived prescriptions of using TIR/UV and FUV-NUV as 
UV dust attenuation indicators by employing a new empirical method. Specifically, these 
calibrations were obtained by fitting a TIR/UV versus FUV-NUV relation, which was de- 
fined by assuming that the UV attenuation estimated from TIR/UV ratio matches that from 
FUV-NUV color, to a nearby star-forming galaxy sample. We have demonstrated that the 



11 When FUV is available, there is no reason to use NUV as a SFR measure. So the FUV-NUV corrected 
NUV luminosity is not examined. 
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combinations of the observed FUV and NUV luminosities with the TIR luminosities (i.e., 
the TIR/UV corrected UV luminosities) are tightly correlated with the Balmer-corrected Ha 
luminosities, and hence can serve as good SFR measures. By referencing to the TIR/UV 
corrected UV luminosities, we used the 25/zm (or 24/im) infrared or 1.4 GHz radio continuum 
luminosities as surrogate of the TIR luminosity to correct for the dust attenuation in the 
UV and found that they are slightly less accurate than the TIR luminosity. The combina- 
tions adopt the form of L(UV) corr = L(UV) obs + a ■ L(IR or radio), and the corresponding 
attenuation in the UV has the form of A uv = 2.51og[l + a ' L ^° r , radio) ], where a is the ap- 

V jobs 

proximation of the product of the inverse of the bolometric correction r/uv and the factor 
^7^-^ (see Section 3). The coefficient a and the dispersions of the combinations relative 
to the Balmer-corrected Ha luminosities are listed in Table 3. 

As we emphasized in Paper I, these composite indices have been calibrated over lim- 
ited ranges of galaxy properties and physical environments. It is therefore important to 
understand the range of observations over which these calibrations are derived, and any 
systematics which were built into our calibrations. It is also useful to compare our results 
with the calibrations presented in the literature by different groups. We also reassess the 
consistency between the SFR derived from our dust attenuation corrected FUV luminosity 
and that based on Balmer-corrected Ha luminosity. 



6.1. Range of Applicability 

In Paper I we have already discussed the ranges of the galaxy properties (e.g., the 
attenuation in Ha, TIR to Ha flux ratio, D n (4000A) and TIR luminosity) covered by our 
samples. Briefly our sample is composed of normal star-forming galaxies in the local universe. 
Neither early-type red galaxies nor dusty starburst galaxies are included in this sample, which 
is further demonstrated by the UV properties presented in this paper. The ranges of IRX 
and FUV-NUV color spanned by our sample galaxies are -0.19 to 2.39 dex and 0.07 to 1.06 
mag respectively. This coverage in IRX and FUV-NUV color introduces uncertainties in the 
best fitted coefficients, as presented in Section 3, due to the lack of constraints from objects 
redder than 1.06 mag and bluer than 0.07 mag. As a result caution should be exercised when 
applying our calibrations to objects outside the above ranges (see Section 6.6). 
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6.2. Comparisons With Calibrations In The Literature 

There exist several calibrations for normal star-forming galaxies in the literature. These 
calibrations were obtained under different assumptions about detailed star formation his- 
tories, dust attenuation curves, and the parameterization of the A FU v ~ IRK relation (see 
below). It is thus necessary to test the consistency between them. We compare our calibra- 
tion for IRX-based FUV attenuation estimates with those commonly used in the literature 
in Figure 13. Here the FUV attenuation (top panel), the residual of the FUV attenuation 
estimated using the calibrations given by other groups relative to that using our calibration 
(middle panel) and the normalized (percentage) residual (bottom panel) are plotted as a 
function of IRX. The black solid line denotes our calibration. The calibrations presented by 
other authors are shown with different line types as denoted in the top panel. Specifically, 
the commonly used calibrations for star- forming galaxies by Buat et al. (2005), Burgarella 
et al. (2005) and Kong et al. (2004) are plotted. We also show the well-known prescription 
for starburst galaxies by Meurer et al. (1999). The shaded region denotes the uncertainty in 
our calibration that is estimated using error propagation of eq. (15), without considering the 
measurement errors in FUV and TIR luminosities. The range of this region in IRX indicates 
the coverage by our sample objects, -0.19 to 2.39 dex (see Section 6.1). Outside this range 
the relation shown by the black solid line is simply the extrapolation of that defined by our 
data. 

Before comparing the calibrations in the literature to our Apuv ~~ IRX relation, it is 
interesting to compare the published relations to each other first. From Figure 13, it can be 
clearly seen that over the whole range of IRX probed by our sample galaxies, the relations 
given by Buat et al. and Burgarella et al. are almost indistinguishable from each other - 
the difference is 0.1 mag at maximum. When compared to the Kong et al. relation, which is 
also derived for normal star-forming galaxies, at low IRX (L(TIR)< 10L(FUV) o b s ), the three 
relations for star-forming galaxies stay close to each other - the difference is no more than 0.1 
mag in the estimated attenuations. But at higher IRX (L(TIR)> 100L(FUV) obs ) the Kong 
et al. relation shows significant difference from the other two, 0.3-0.5 mag. Interestingly, 
the difference between the starbursts relation and the Buat et al. and Burgarella et al. 
relations is not dramatic (no more than 0.5 mag), given that they are derived for different 
types of galaxies. This comparison tells us for galaxies with —0.19 < IRX < 2.39, when the 
attenuation prescriptions were used blindly, regardless of galaxy types, the uncertainty in 
the UV attenuation would be less than 0.5 mag. 

Since our sample galaxies are mostly normal star-forming galaxies in the local universe 
(see Section 6.1), it is expected that our calibration is more similar to those calibrated 
for star-forming galaxies than for starbursts. Figure 13 confirms this expectation. The 
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starburst curve obtained by Meurer et al. systematically overestimates the Apuv compared 
to our calibration, although the difference is not dramatic, still within la uncertainty of our 
calibration. The difference between starbursts and star-forming galaxies is expected because 
the mean age of the dust heating population should be systematically lower in starburst 
galaxies than in normal star-forming galaxies. Overall our calibration is very close to those 
for star-forming galaxies derived by Burgarella et al. and Buat et al. over the whole range 
of IRX spanned by our sample. It is worth noting that the prescriptions given by Burgarella 
et al. and Buat et al. are in polynomial forms, different from the form presented in eq. (2), 
which contributes to part of the differences. Among the three star-forming relations that we 
explored here, that of Kong et al. (2004) appears to deviate the most, as seen in Figure 13. 
At < IRX < 1.8 the consistency between Kong et al. relation and ours is good, but at 
IRX < and IRX > 1.8 the relation in Kong et al. over- and under-estimate A FU v at > la 
level, respectively. 

We compare our calibration on Apuv versus FUV-NUV relation with those in the lit- 
erature in Figure 14. Similar to Figure 13, the Apuv versus FUV-NUV relation is plotted 
in the top panel, the difference in Apuv between the prescriptions in the literature and ours 
is shown in the middle panel, and the normalized (percentage) residual is shown in the bot- 
tom panel. The shaded region denotes the uncertainty in our calibration over the range of 
FUV-NUV color covered by our sample galaxies. As before the relation obtained in Section 
3 is plotted with a black solid line and the calibrations presented by other authors are rep- 
resented by different types of lines as denoted in the top panel. In this case, the relations 
plotted include the starburst galaxy calibration by Meurer et al. (1999), the calibrations for 
spectroscopically and NUV-r color selected star-forming galaxies by Salim et al. (2007), the 
relation by Treyer et al. (2007) for star-forming galaxies selected spectroscopically and the 
relation by Seibert et al. (2005) for a wide assortment of galaxy types. 

It is obvious that at a given FUV-NUV color, the starburst line derived by Meurer et 
al. (1999) overestimates the Apuv of normal star-forming galaxies dramatically compared 
to others. This is consistent with the well-known deviation of star-forming galaxies from 
starbursts in IRX-/3 diagram (e.g., Kong et al. 2004; Cortese et al. 2006; Dale et al. 2007; 
Johnson et al. 2007b). The relations defined by star- forming galaxies stay close to each other 
but with different slopes and intercepts. Most of the differences shows up at FUV-NUV < 1. 
For FUV-NUV > 1, the relations converge except for those by Salim et al., which have 
constant A F uv for galaxies redder than 0.90 and 0.95 in FUV-NUV color for NUV-r and 
spectroscopically selected galaxies respectively. As discussed in Section 4.3, the calibration 
derived by Treyer et al. (2007) underestimates the attenuation by 0.4 mag at FUV-NUV 
~ mag with the discrepancy decreasing as the object gets redder in FUV-NUV color. At 
FUV-NUV=1.7 mag, the Treyer et al. calibration gives consistent result with ours. The 
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calibration by Seibert et al. (2005) is the most consistent one with ours, and the estimated 
Afuv from their relation agrees with that from ours within ±la for the range of FUV-NUV 
0.2-1.7 mag. For bluer color (0 - 0.2 mag), the difference is a bit larger but always smaller 
than 0.2 mag in absolute scale. In terms of relative difference (bottom panel), both Treyer 
et al. and Seibert et al. relations are consistent with ours within ±1(7 at FUV-NUV > 0.5, 
with increasing discrepancy at bluer color. For Salim et al. calibrations, both the absolute 
and the relative differences are large, far beyond ±1<t, for galaxies redder than ~ 0.90 in 
FUV-NUV because of the plateau in A FUV defined in these relations, which is hard to be 
tested using our sample because of the limited range in FUV-NUV color covered by our 
galaxies (as indicated by the shaded region). 



6.3. Comparison with Dust Attenuation Laws 

As shown in Section 3, the scale parameter spuv i n equation (9) represents the slope 
of the UV part of the attenuation curve. We derive /cfuvAnuv from spuv and compare it 
with those derived from the commonly used dust attenuation laws published by Calzetti et 
al. (2000; originally from Calzetti et al. 1994;) and Chariot & Fall (2000). The attenuation 
curve by Calzetti et al. (1994) was obtained directly from the FUV-to-near-IR spectra 
of a sample of UV-selected starburst galaxies. Chariot & Fall (2000) used a similar but 
larger sample to derive their effective attenuation curve by tuning their model parameters 
to account for the distribution of the galaxies in the IRX, Ha/H/? ratio and EW(Ha;) versus 
(3 diagrams. The estimated values of /cfuvAnuv from our calibration, Calzetti et al. (2000) 
and Chariot & Fall (2000) are 1.35±0.06, 1.24 and 1.32, respectively, which means that in 
the UV part of the dust attenuation curve, the slope of the attenuation curve defined by our 
sample is consistent with that derived using Chariot & Fall (2000) attenuation law within 
the uncertainties but marginally steeper (at ~ 2a level) than that estimated from Calzetti et 
al. (1994;2000). The slope of our attenuation curve also suggests that the global FUV and 
NUV attenuations show no evidence for the effects of a strong 2175A bump (e.g., O'Donnell 
1994). For example, the O'Donnell Milky Way extinction curve has a /cfuv/^nuv ratio close 
to 0.9, which is much smaller than 1.35±0.06 defined by our galaxies. 

Based on the attenuation in the FUV estimated using our prescriptions, we can study 
the Apuv ~~ Ahq relation. Figure 15 shows the IRX-derived (left panel) and FUV-NUV 
derived (right panel) Apuv as a function of Ha/H/3-derived Ah q . The solid line in this 
figure is derived from eqs. (14) and (16), which present the relations between FUV-NUV 
color and Ah q and Afuv respectively. For comparison, the relation derived from Calzetti 
(2001) law after adapted to GALEX FUV waveband, which has the form A FUV = 1.82A Ha , 
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is plotted as the dotted line. As can be seen from the figure, at a given Ha attenuation, 
Calzetti's law gives a lower FUV attenuation estimate than ours and the difference becomes 
more significant for galaxies that are more attenuated. Quantitatively speaking, however, 
the slope of the solid line, which determined from our data, is 2.06±0.28, consistent with 
that given by Calzetti's law within the error. 

6.4. Comparisons Made For NUV 

We performed similar comparisons for NUV with the calibrations published in the lit- 
erature. When compared to the calibrations in the literature, similar to the results we saw 
for FUV, our calibration is in good agreement with the one given by Burgarella et al. (2005) 
within 0.1 mag for the whole range of Log[(TIR)/L(NUV)]. The Anuv/Aho ratio defined by 
our data is 1.52±0.56, coincident with the predicted value from Calzetti (2001)'s law within 
the error. 



6.5. Dust Heating From Old Stars 

It is now well-known that a significant fraction of IR emission comes from dust heating 
by old stellar populations (e.g., Sauvage & Thuan 1992; Bell 2003 and references therein). In 
Paper I, by comparing our calibration for Ha+TIR as a SFR indicator with stellar population 
synthesis models, we found that for our sample objects up to 50% of the TIR emission could 
be from stars older than 100 Myr. Are the results presented here consistent with those 
estimates? 

As mentioned in Section 3, the coefficient apuv i n ec L- (2) is the product of the inverse 
of the bolometric correction ?7fuv an d the factor ^-j^pFp- In Paper I, we have empirically 
derived f by combining the TIR luminosity and the multi-wavelength SEDs of a subset of 
MK06 galaxies with GALEX, SDSS and 2MASS observations. Briefly, eq. (6) was used to 
estimate f from the TIR and the bolometric luminosities, which were derived by adding the 
TIR luminosity to the integrated luminosity from the FUV to K band. After substituting 
the values of tfuv and f into riy-^, we find that the median of ^r^pjy^ is around 1.5. 
Then the bolometric correction 1/t/fuv is estimated to be ~ 3.3 by comparing ofuv with 

(1-C~ T FUV) 
(1-0--) • 

In order to estimate the effect of dust heating from old stars, we compare the bolometric 
correction derived from our data to those from stellar population synthesis models. We found 
that the bolometric correction defined by our data is more than two times greater than the 



-26 - 



model prediction for a 100 Myr old stellar population with constant star formation history. 
This implies that stars older than 100 Myr must contribute to the dust heating. The spectral 
synthesis modelling using STARBURST99 indicates that the bolometric correction for a 10 
Gyr old stellar population with constant star formation history, solar metallicity and either 
a Salpeter or a Kroupa IMF is only 15%-20% lower than the estimated value above. In 
other words, the mean property of our sample objects approximates a stellar population 
that has been constantly forming stars for 10 Gyr, which is exactly what we inferred from 
the analyses of Ha+TIR in Paper I. It is worth noting that the consistent inference drawn 
from Paper I and this paper ensure us the robustness of our results. However, this estimate 
of the average stellar population of our sample galaxies conflicts with the dust-free FUV- 
NUV color derived in Section 3, which is ~ 0.04 mag redder than a 10 Gyr old galaxy with 
constant star formation history. All the estimates made here are very rough and they should 
only be used as a guide to our understanding of the underlying physics shown by the data 
qualitatively. More careful modelling is needed for a quantitative understanding. 



6.6. Systematic Uncertainties 

As discussed in Section 6.1, our calibrations were derived over limited ranges of galaxy 
properties. So it is important to understand possible systematics which were built into our 
prescriptions. 

First of all, prescriptions calibrated for different types of galaxies are different. This is 
not surprising because different types of galaxies show different IRX-/3 (i.e., FUV-NUV color) 
relations. As shown in Section 6.2, the IRX-based and the FUV-NUV color-based attenuation 
correction prescriptions for star-forming galaxies differ from those for starburst galaxies. 
Another two extreme examples are dwarf irregulars and (Ultra) Luminous Infrared Galaxies 
((U)LIRGs). Dale et al. (2009) examined the location of the Local Volume Legacy (LVL) 
sample galaxies in the IRX-/3 diagram. Nearly two-thirds of LVL galaxies are dwarf /irregular 
systems. Those authors showed that the majority of LVL galaxies fall below the IRX-/3 
relation as defined by normal star-forming galaxies, and they cluster in an area with blue 
FUV-NUV color and low IRX (see also Lee et al. 2009). The range in IRX spanned by 
the LVL galaxies, for a given FUV-NUV color, is large - more than an order of magnitude. 
For (U)LIRGs, Howell et al. (2010; see also Goldader et al. 2002; Calzetti 2001) showed 
that (U)LIRGs in the Great Observatories ALL-sky LIRG Survey (GOALS) lie above the 
starburst relation, and the correlation between IRX and FUV-NUV color is not strong either. 
Similar to dwarf/irregulars, at a given FUV-NUV color, the (U)LIRGs span more than an 
order of magnitude in IRX. One assumption we made in our calibration method is that 
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the attenuation derived from IRX matches that estimated from FUV-NUV color, which is 
reasonable for our normal star-forming galaxy sample given the relatively tight correlation 
between IRX and FUV-NUV color, as shown in Figure 3. However, it breaks down for dwarf 
irregulars and (Ultra) Luminous Infrared Galaxies given the large range of IRX spanned by 
these galaxies at a given FUV-NUV color. Therefore no attempt was made to calibrate an 
IRX-based and FUV-NUV color-based attenuation correction prescriptions for these galaxies. 

Apart from different galaxy samples, the differences between different calibrations can be 
traced to the different methodologies. Previous calibrations have been based either on pure 
model assumptions, or on a comparison of observational data with models. The difference 
in the form of the parameterization of the prescriptions also contributes to the differences. 
As we can see from the comparisons in Section 6.2, the differences in methodology and 
parameterization (Buat et al. 2005; Burgarella et al. 2005 and this work) only account for a 
small portion of the differences, which cannot even be distinguished within ±la uncertainty, 
and the main contributors to the different relations are the star formation history and the 
dust attenuation curve that includes the effect of the dust / star distribution, as suggested by 
other authors (e.g., Chariot & Fall 2000; Kong et al. 2004; Burgarella et al. 2005). 

It is difficult to test the impact of dust attenuation curves on attenuation calibrations. 
However, it is possible to estimate the effect of dust-heating stellar populations. As dis- 
cussed in Section 6.5, for our sample galaxies, up to half of the TIR emission could be from 
dust heated by stars older than 100 Myr. By contrast, in starburst galaxies, the dust is 
dominantly heated by stellar populations younger than 100 Myr (Meurer et al. 1999). This 
difference in dust-heating stellar populations of starbursts and our star-forming galaxies 
results in the difference in attenuation estimates, 0.3 mag at maximum over —0.19 < IRX 
< 2.39, as shown in Figure 13. The effects on the UV color-based corrections are much larger 
than on IRX-based prescriptions. The IRX-based prescriptions are surprisingly consistent 
with each other. However different FUV-NUV based attenuation schemes have by compari- 
son enormous inconsistencies. This underscores our concerns about color-based attenuation 
corrections based on scatter in our own data. 

In principle, by comparing modelled SEDs with multi-wavelength data, the differences 
in star formation histories and dust attenuation curves can be taken into account on a 
galaxy-to-galaxy basis. However, the dust attenuations derived this way and the resulting 
relations of Apuv versus IRX or Apuv versus FUV-NUV color are dependent on the adopted 
stellar population synthesis models and the dust attenuation curves. On the other hand, the 
requirement of the availability of large libraries of models and multi-wavelength data makes 
this method less practical in many cases. 

By contrast, our approach is more straightforward and less mo del- dependent. Our cal- 
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ibrations rely almost entirely on the data themselves and hence can serve as independent 
checks of other model-dependent methods. Although we assumed an average stellar popula- 
tion and dust attenuation curve (as discussed in Section 3), they do not introduce systematics 
into the calibrations and are valid on a statistical basis. Therefore, our calibrations should 
be robust. They also tend to be consistent with most of the published calibrations for normal 
star-forming galaxies. 

The last possible systematic error comes from the TIR luminosities we adopted in our 
calibrations. We used TIR luminosities based on IRAS observations, which are systematically 
lower than those based on Spitzer/MIPS measurements by 24% ± 23% on average (Section 
2 and Figure 1 in Paper I). If we simply substitute the IRAS TIR luminosities by the MIPS 
TIR luminosities, the coefficient in front of IRX in eq. (15) will be lower by 24%, i.e., 
changes to 0.37, which is within the la uncertainty of the coefficient. So if one uses MIPS 
TIR luminosities to estimate the FUV attenuation using our IRX-based prescription (i.e., 
eq. (15)), the uncertainty will be less than 0.3 mag. 

In summary, our IRX-based attenuation correction has an uncertainty of < 0.3 mag, 
when it is applied to starburst galaxies or to galaxies with MIPS based TIR luminosities. 
By contrast, if one uses our FUV-NUV color prescription to estimate the FUV attenuation 
for starbursts, he can underestimate the FUV attenuation by more than 1-2 mag. But for 
dwarf irregulars or (U)LIRGs, our prescriptions may suffer from severe problems. 

6.7. The SFR Matching Method 

In the introduction, we mentioned that Treyer et al. (2007) derived an A FUV - FUV- 
NUV relation by matching the FUV-based and Ha-based SFRs. The attenuation derived 
this way is dependent on the intrinsic ratio of the FUV to Ha luminosity, which in turn 
is dependent on model assumptions. Here we use a similar method to re-calibrate the 
A.fuv~IRX relation in order to evaluate the model dependence of the calibration obtained 
this way. We converted the Balmer-corrected Ha luminosity to a Balmer-corrected FUV 
luminosity using a theoretical FUV/Ha ratio and then force the combined observed (i.e., 
without internal dust attenuation correction) FUV and TIR luminosity to match the Balmer- 
corrected FUV luminosity so that the coefficient a F uv i n eq. (2) can be obtained. We used 
the same set of models as those used in Figure 5 and Figure 8. The model predictions and the 
best-fit coefficients a F uv are listed in Table 2 (For reference the corresponding quantities for 
NUV are also listed.). Consistent with what we saw from Figure 5, Table 2 shows that our 
calibration derived in Section 3 is in excellent agreement with that predicted by assuming a 
constant star formation history with Kroupa IMF and age 100 Myr. 
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The FUV/Ha ratio changes at most by 10% when only the effect of a single parameter 
is taken into account, and accordingly the resulting coefficient ofuv varies by 11%. The 
change in the coefficient directly affects the estimated attenuation from IRX. This impact 
is shown in Figure 16. It is clear that assuming a Salpeter or a Kroupa IMF does not affect 
the calibrations dramatically - at maximum ~ 0.07 mag over the range of IRX=0 - 3. The 
calibrations are also affected by assuming either a 100 Myr old or a 1 Gyr (or older) old stellar 
population - at maximum ~ 0.11 mag over the range tested. However the most widely used 
K98's prescriptions give larger discrepancy - -0.2 mag at maximum when IRX=3. This large 
discrepancy with K98's prescriptions is probably due to the use of the old stellar evolutionary 
models, which leads to an over-estimate of Ha luminosity by 11% compared to that from 
STARBURST99 under the same set of assumptions of IMF, metallicity and star formation 
history. 

From these tests, we can draw the following conclusions: (1) The calibration established 
using our new method is consistent with that based on SFR matching method for a constant 
star formation history with Kroupa IMF and age 100 Myr. (2) Compared to our prescription, 
the adoption of a Salpeter IMF and/or an older age both leads to an overestimate in Apuv 
by up to 0.1 mag over the range of IRX=0 - 3. (3) Relative to our calibration, the widely 
used K98's SFR prescriptions can underestimate Apuv by up to 0.2 mag (at IRX=3). 

7. SUMMARY 

We have calibrated log [L(TIR)/L(FUV) obs ] (i.e., IRX) and the FUV-NUV color as 
FUV dust attenuation indicators using a nearby star-forming galaxy sample studied by 
Moustakas & Kennicutt (2006). The combinations of 25/im and 1.4GHz radio continuum 
luminosities with the observed FUV luminosities were also empirically calibrated to probe 
the dust attenuation corrected FUV luminosities. Similar calibrations were derived for NUV 
band as well. The coefficients in these calibrations are summarized in Table 3 (see Section 
6). These prescriptions provided in Table 3 can be used to derive attenuation-corrected star 
formation rates. Our main results can be summarized as follows. 

1. The IRX-corrected FUV luminosities based on our new calibration, i.e., 

L(FUV) corr = L(FUV) obs + (0.46 ± 0.12)L(TIR), 

show tight and linear correlation with the Balmer decrement corrected Ha luminosities, with 
a rms scatter of ±0.09 dex (±0.23 mag). Statistically speaking, their ratios are consistent 
with the model predictions from STARBURST99 by assuming a constant star formation 
rate over 100 Myr and solar metallicity. This consistency applies whether a Salpeter or a 
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Kroupa IMF with lower and upper mass limits of 0.1 and 100 M is adopted. Furthermore, 
these dust attenuation corrected FUV/Ha ratios do not show any trend with other galactic 
properties over the ranges covered by our sample objects. These results suggest that lin- 
ear combinations of TIR luminosities and the observed FUV luminosities (without internal 
attenuation corrections) are excellent star formation rate tracers. 

2. The FUV-NUV corrected FUV luminosities are broadly correlated with the Balmer- 
corrected Ha luminosities. But they do not trace each other linearly and the rms scatter 
is large - ±0.23 dex (±0.58 mag), which is ~ 2.5 times larger than the case for IRX- 
corrected FUV luminosities. In addition, the attenuation corrected FUV/Ha ratios show 
correlations with a few other galactic properties, which will introduce systematic errors into 
the attenuation estimates in different galaxy samples. We confirm others' findings that FUV- 
NUV color is not a good dust attenuation indicator for normal star forming galaxies though 
it can be used with caution when IR data are not available. 

3. Linear combinations of 25/im and 1.4GHz radio continuum luminosities with the 
observed FUV luminosities can be used as surrogates of IRX-corrected FUV luminosities to 
trace the attenuation-corrected star formation rates. Their correlations with the Balmer- 
corrected Ha luminosities are slightly less tight than the IRX-corrected FUV luminosities - 
0.13 dex and 0.14 dex, and a non-linearity is shown in the correlation between FUV±1.4GHz 
and Balmer-corrected Ha luminosities with a slope of 1.09 ± 0.02. 

4. Our calibrations are tested on normal star-forming galaxies with Balmer-based A Hq 
= 0-2.5 mag, IRX = -0.19 - 2.39 and FUV-NUV = 0.07 - 1.06 mag. The A FUV - IRX 
relation shows overall consistency with prescriptions for star-forming galaxies by Burgarella 
et al. (2005) and Buat et al. (2005) within ±ltr uncertainty. While our A FUV - FUV-NUV 
color relation is in good agreement with that presented by Seibert et al. (2005). When 
our IRX-based attenuation correction is applied to starburst galaxies or to galaxies with 
MIPS based TIR luminosities, it has an uncertainty of < 0.3 mag. By contrast, if our 
FUV-NUV color prescription is used to estimate the FUV attenuation for starbursts, it can 
underestimate the FUV attenuation by more than 1-2 mag. 

5. The estimated A F uv using our prescriptions is related to the Balmer-based A Ha in a 
way consistent with that given by Calzetti (2001) 's law within the error. 

6. Our analyses show that the combination of the NUV and the TIR luminosities can 
serve as a rather good SFR indicator for our sample galaxies, but may become problematic 
in galaxies with non-constant star formation histories in the past ~ 1 Gyr. 
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Fig. 1. — Comparison of photometry from GALEX AIS images with that from images with 
longer exposure time for MK06 sample at FUV (left) and NUV (right) wavebands. 
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A Ha ( ma §) 

Fig. 2. — FUV-NUV color as a function of attenuation in Ha line derived from the Balmer 
decrement ratio for MK06 sample (solid circles) and SINGS sample (open circles), with the 
corresponding UV continuum slope /3 labelled on the right axis. The solid line is the bisector 
fitting to the MK06 sample. The error bars in the bottom-right corner denote median errors 
for the MK06 sample (right) and the SINGS sample (left). 
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Fig. 3.— IRX (i.e. log [L(TIR)/L(FUV)]) versus FUV-NUV color for the MK06 sample 
(solid circles) and the SINGS sample (open circles), with j3 labelled on the top axis. The 
open squares with error bars represent the median, lower (25%) and upper (75%) quartiles 
in bins of width 0.2 mag in FUV-NUV color. The solid line, which has the form defined 
in eq. (13), is the best fit to the MK06 sample. The error bars in the bottom-right corner 
denote median errors for the MK06 sample (right) and the SINGS sample (left). 
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Log [L(Ha) corr ] (ergs s" 1 ) 

Fig. 4. — Observed FUV luminosities without correction for internal dust attenuation versus 
Balmer decrement ratio corrected Ha luminosities for MK06 sample (solid circles) and SINGS 
sample (open circles). The color-coded lines represent the predicted relations (in the absence 
of dust attenuation in the UV) by different SFR prescriptions, same as Figure 5. The error 
bars in the bottom-right corner denote median errors for the MK06 sample (right) and the 
SINGS sample (left). 
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Fig. 5. — left panel: IRX-corrected FUV luminosities as a function of Balmer decrement ratio 
corrected Ha luminosities for MK06 sample (solid circles) and SINGS sample (open circles). 
The lines overplotted represent the predicted relations by different SFR prescriptions (see 
text), which are the K98 prescription (red short-long dashed line) and models constructed 
using STARBURST99 by assuming a constant star formation history, solar metallicity for 
Kroupa IMF at 100 Myr (blue short dashed line), Kroupa IMF at 1 Gyr (blue long dashed 
line), Salpeter IMF at 100 Myr (magenta short dashed line) and Salpeter IMF at 1 Gyr 
(magenta long dashed line). The error bars in the bottom- right corner denote median errors 
for the MK06 sample (right) and the SINGS sample (left), right panel: Histograms of the 
ratios of the IRX-corrected FUV luminosities to Balmer-corrected Ha luminosities for MK06 
sample (solid line) and SINGS sample (dotted line) in log scale. The vertical lines represent 
the predicted values by different SFR prescriptions, as shown in left panel. This panel is 
ffiven to show the differences between the models more clearlv 
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Fig. 6. — The logarithm residuals of the attenuation-corrected FUV luminosities from 
IRX relative to the Balmer-attenuation-corrected Ha luminosities as functions of Balmer- 
attenuation-corrected Ha luminosities (top-left panel), Balmer-attenuation-corrected Ha lu- 
minosities per unit area (top-right panel), attenuation in Ha calculated from Ha/H/3 ratio 
(bottom- left) and axial ratio (inclination) b/a (bottom- right). The label at the top of the 
top-left panel shows the corresponding SFR according to the L(Ha) - SFR relation given by 
K98. The solid circles represent MK06 galaxies while the open circles denote SINGS galax- 
ies. The dashed line denotes the predicted value by the STARBURST99 synthesis model for 
a constant star formation history, solar metallicity and Kroupa IMF at 100 Myr. 
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Fig. 7. — Similar to Figure 6 but with residuals plotted as functions of integrated Equivalent 
Width at Ha (top- left panel), 4000 A break (D n (4000A)) (top-right panel), gas-phase oxy- 
gen abundance (12+log(0/H)) (bottom-left panel) and FIR color (/ iy (60/im)// iy (100/im)) 
(bottom- right panel). See Figure 6 for the explanation of the symbols. 
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Fig. 8. — left panel: FUV-NUV color-corrected FUV luminosities as a function of Balmer 
decrement ratio corrected Ha luminosities for MK06 sample (solid circles) and SINGS sam- 
ple (open circles). The lines overplotted represent the predicted relations by different SFR 
prescriptions, same as those in Figure 5. The error bars in the bottom-right corner denote 
median errors for the MK06 sample (right) and the SINGS sample (left), right panel: His- 
tograms of the ratios of the FUV-NUV color-corrected FUV luminosities to Balmer-corrected 
Ha luminosities for MK06 sample (solid line) and SINGS sample (dotted line) in log scale. 
The vertical lines represent the predicted values by different SFR prescriptions, as shown in 
left panel. This panel is to show the differences between the models more clearly. 



-44 - 



Log [SFR(Ha) corr ] (M yr"«) 



-2-1012 




39 40 41 42 43 37 38 39 40 41 

Log [L(Ha) corr ] (ergs s" 1 ) Log L(Ha) oorr /area (ergs s"> kpc" 2 ) 




Fig. 9. — The logarithm residuals of the attenuation-corrected FUV luminosities from 
FUV-NUV color relative to the Balmer-attenuation-corrected Ha luminosities as func- 
tions of Balmer-attenuation-corrected Ha luminosities (top-left panel), Balmer-attenuation- 
corrected Ha luminosities per unit area (top-right panel), attenuation in Ha calculated from 
Ha/H/9 ratio (bottom-left) and axial ratio (inclination) b/a (bottom-right). The label at 
the top of the top-left panel shows the corresponding SFR according to the L(Ha) - SFR 
relation given by K98. The solid circles represent MK06 galaxies while the open circles de- 
note SINGS galaxies. The dashed line denotes the predicted value by the STARBURST99 
synthesis model for a constant star formation history, solar metallicity and Kroupa IMF at 
100 Myr. 
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Fig. 10. — Similar to Figure 9 but with residuals plotted as functions of integrated Equiv- 
alent Width at Ha (top- left panel), 4000 A break (D n (4000A)) (top-right panel), gas-phase 
oxygen abundance (12+log(0/H)) (bottom-left panel) and FIR color (/ i ,(60//m)// i ,(100//m)) 
(bottom-right panel). See Figure 9 for the explanation of the symbols. 




Fig. 11. — Linear combinations of the observed FUV luminosities with 25/im luminosities 
(left panel) and 1.4GHz luminosities (right panel) compared to the Balmer-corrected Ha 
luminosities, with the scaling coefficient derived by matching the combined FUV and 25fim 
or 1.4GHz luminosities with the IRX-corrected FUV luminosities. The solid circles represent 
MK06 galaxies, the open circles denote SINGS galaxies with the integrated measurements 
while the open squares denote the central 20" x 20" regions of SINGS galaxies. The dashed 
line denotes the predicted value by the STARBURST99 synthesis model for a constant star 
formation history, solar metallicity and Kroupa IMF at age 100 Myr. The error bars in 
the bottom-right corner denote median errors for the MK06 sample (right) and the SINGS 
sample (left). 
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Fig. 12. — left panel: Combined NUV and TIR luminosities as a function of Balmer decre- 
ment ratio corrected Ha luminosities for MK06 sample (solid circles) and SINGS sample 
(open circles). The error bars in the bottom- right corner denote median errors for the MK06 
sample (right) and the SINGS sample (left), right panel: Histograms of the ratios of the 
TIR/NUV-corrected NUV luminosities to Balmer-corrected Ha luminosities for MK06 sam- 
ple (solid line) and SINGS sample (dotted line) in log scale. The color-coded lines represent 
the same model predictions as those in Figure 5. 
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Fig. 13. — Comparison of the Afuv vs. IRX (i.e., log [(L(TIR)/L(FUV)] relation derived in 
this paper to those presented by others. The comparison is made in three different ways: 
the A FU v itself (top), the difference in A FU v relative to our calibration (middle) and the 
normalized difference (bottom) are plotted as a function of IRX. The solid line denotes our 
result. The shaded region indicates the uncertainty in our calibration over the range of IRX 
covered by our sample galaxies. 
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Fig. 14. — Top: Comparison of the Apuv versus FUV-NUV relation obtained in this work to 
those presented by others. Middle: The difference in Apuv between the relations given in the 
literature and that in this paper versus FUV-NUV. Bottom: The normalized residual versus 
FUV-NUV. The solid line denotes our result. The shaded region denotes the uncertainty in 
our calibration over the range of FUV-NUV color spanned by our sample galaxies. 
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Fig. 15. — IRX (left panel) and FUV-NUV (right panel) derived dust attenuation in FUV 
as a function of the dust attenuation in Ha line. The solid line represents the relation 
defined by our sample galaxies and the dotted line denotes that given by Calzetti (2001), 
with the difference in E(B-V)* and E(B-V) gas taken into account and adapted to GALEX 
FUV waveband. The error bars in the bottom-right corner denote median errors for the 
MK06 sample (right) and the SINGS sample (left). 
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Fig. 16. — Top: Apuv versus IRX (i.e., log [L(TIR)/L(FUV)]) relation with coefficient a 
derived using the intrinsic FUV to Ha luminosity ratio based on different stellar population 
assumptions, i.e., IMFs and ages (see text), a is defined by A FUV = 2.5 log(l + ofuv ■ 10 IRX ) 
and is the approximation of the product of the inverse of the bolometric correction //fuv and 
the factor ■ Bottom: The difference in A FUV for different model assumptions. 
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Table 1. GALEX UV photometry for MK06 sample 



Name 


Distance 


FUV 


NUV 




(Mpc) 


(mjy) 


(mjy) 


ARP256 
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q qq+n 1 1 
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00 41 _Lf| «0 
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J- . (jiJ-l—VJ . -L<J 
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1 1 oo-Un ^9 


1 q c;4-un ^4 


TCI 623 


84.1 


k nO+0 28 


7 01 _Ln 99 


iunn 03 04 ni 4 


139.9 




f| c:o_Lf| f)i 


ncco6Q5 


135.4 




1 40+0 04 

± . lUIU.U^ 


NCC0877 


54.6 


5.11±0.47 


7 e;fi_Lf) 4f| 

1 - r M 1 1 11. '-H ) 




43.4 


1401 _|_n 

14.Z1IU.UU 




ncco958 


79.5 


l.vU _L VJ . KJ 1 


9 4fi+n 07 


NCC0959 


10.9 


A c:O + 21 


a ki _i_n 1 8 




59.8 


1 07-Un 98 
i.y 1 icu.ZrO 


q 90-kO 9^ 


1 1 VJv J- woo 


9.8 


c s«+0 49 


q 01 _i_n Qq 


noci ns4 


20.0 


1 q 77+n fi4 
±0.1 f iu.ui 


99 c;9-un fi9 


IN V_J V_> J_ IrrU 


21.2 


n 7fi-un 4^ 


1 1 7fi-Un 33 


Nnci 1 56 


7.8 


24.71dzl.14 


qq 4q+n S4 


i> v.j v. iojy 


28.0 


J. X . OIIU.UO 


1 q 7Q-Ufi 49 


1 1 VJVj J-tJOtJ 


21.2 


1 7 qo+O 83 


94 7fi+n fiQ 


TR AS0335Q+1 523 


147.8 


n 71 +n 36 

U. 1 J — l_U . vj 


1 25+0 27 




29.7 


y.uyiu.uu 


1 4 QQ-un ^4 


TTCC02Q82 


74.5 


n fiQ-U9 7S 


n Qf;+i on 

VJ . iS'J _l_ l.l/U 


Nnci 569 


2.0 


379 4fi-M 7 1fi 


«71 91 _L1 Q R9 


NCC21 39 


26.7 


1 q QQ-un 79 


9fl 9Q-un (SI 


NCC2337 


7.9 


4 «4+n 31 


fi 97+0 9^ 


NCC941 5 

IN Vjt vj^l J_ tj 


57.2 




q 41 _Ln 9Q 


NCC2500 

1 N vJ V— AtJVJVJ 


11.0 


1 q qq+o (S9 


1 k 1 7+n 49 


IN \j z,oo 1 


6.9 


q D7+0 49 


1 pi 77-un 30 

1U. ( 1 IU.OU 




OD. / 


1.0z±U.Uo 


Z.40±U.U / 


NGC2903 


8.9 


34.67±1.60 


60.86+1.68 


NGC3239 


10.5 


14.91±0.69 


26.41+0.73 


NGC3265 


25.0 


0.59±0.03 


0.97+0.03 


UGC05720 


26.9 


4.05±0.19 


5.05+0.14 


NGC3344 


7.9 


35.41±1.63 


49.66+1.37 


NGC3353 


19.1 


4.95±0.23 


6.10+0.17 


NGC3367 


49.3 


10.21±0.52 


15.56+0.47 


ARP270 


30.3 


17.38±0.80 


22.15+0.61 


NGC3432 


10.4 


15.66±0.75 


22.68+0.66 


NGC3442 


31.6 


2.91±0.15 


3.94+0.13 


NGC3521 


8.6 


17.09±0.79 


34.79+0.96 


NGC3600 


13.3 


2.44±0.11 


3.54+0.10 


IC0691 


23.7 


1.17±0.08 


1.58+0.07 


NGC3726 


19.8 


20.25±0.99 


28.26+0.83 


NGC3738 


4.9 


11.49±0.56 


13.85+0.40 


NGC3769 


19.8 


3.98±0.23 


5.89+0.21 
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Table 1 


—Continued 




Name 


Distance 


FUV 


NUV 




(MpcJ 


(mjy) 


(mjyj 


UGC06665 


87.4 


2.46±0.11 


3.41±0.10 


NGC3870 


19.8 


3.33±0.18 


4.02±0.13 


NGC3893 


15.5 


13.80±0.69 


25.50±0.74 


NGC3928 


16.9 


2.47±0.15 


3.72±0.14 


NGC3949 


19.8 


13.60±0.66 


20.83±0.60 


NGC4004 


55.6 


2.49±0.12 


3.61±0.10 


ARP244 


25.0 


33.06±1.52 


46.56±1.29 


NGC4051 


15.5 


15.22±0.76 


23.94±0.72 


NGC4062 


11.4 


3.96±0.20 


7.13±0.20 


NGC4085 


19.8 


1.40±0.12 


2.70±0.13 


NGC4088 


19.8 


11.37±0.60 


18.43±0.58 


NGC4096 


10.8 


8.04±0.37 


12.58±0.35 


NGC4100 


19.8 


4.24±0.31 


8.24±0.32 


NGC4214 


2.9 


52.74±2.44 


64.17±1.78 


NGC4218 


19.8 


3.26±0.16 


4.05±0.12 


NGC4254 


12.9 


31.13±1.46 


56.30±1.57 


NGC4303 


10.6 


43.55±2.01 


62.09±1.72 


NGC4384 


42.8 


1.99±0.13 


3.06±0.12 


NGC4414 


17.7 


7.46±0.35 


14.08±0.39 


NGC4605 


5.2 


21.53±0.99 


31.48±0.87 


NGC4618 


8.9 


23.61±1.09 


30.54±0.84 


NGC4625 


10.2 


3.88±0.18 


5.25±0.14 


NGC4651 


20.6 


9.52±0.44 


13.98±0.39 


NGC4656 


8.7 


44.27±2.04 


48.63±1.34 


NGC4666 


27.6 


6.14±0.35 


9.15±0.34 


NGC4670 


23.2 


7.82±0.37 


9.38±0.26 


NGC4713 


13.7 


13.05±0.61 


17.91±0.50 


NGC4900 


9.0 


9.86±0.47 


15.42±0.43 


NGC5014 


23.8 


0.86±0.05 


1.55±0.05 


MRK0066 


95.6 


1.36±0.09 


1.60±0.06 


NGC5194 


8.0 


123.02±5.67 


196.70±5.43 


ARP240 


104.9 


4 1 1 _|_n 25 




NGC5253 


3.1 


41.14±1.90 


54.93±1.52 


NGC5430 


49.0 


2.57±0.16 


4.05±0.14 


NGC5607 


112.2 


1.34±0.11 


2.46±0.10 


NGC5591 


117.6 


0.68±0.07 


1.28±0.08 


NGC5653 


58.7 


1.56±0.11 


3.02±0.11 


NGC5676 


37.7 


3.92±0.31 


8.09±0.31 


UGC09560 


23.7 


3.30±0.15 


3.52±0.10 


IC1076 


94.6 


0.88±0.10 


1.54±0.09 


NGC5996 


54.6 


5.63±0.26 


7.93±0.22 


NGC6052 


74.8 


4.99±0.30 


7.60±0.24 


NGC6090 


130.2 


1.15±0.05 


1.99±0.05 


NGC6926 


87.6 


3.08±0.54 


5.34±0.45 


IC5179 


50.3 


3.01±0.18 


5.52±0.20 
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Table 1 — Continued 



Name 


Distance 


FUV 


NUV 




(Mpc) 


(mjy) 


(mjy) 


NGC7316 


79.0 


2.58±0.13 


3.95±0.11 


CGCG453-062 


106.0 


0.13±0.02 


0.35±0.02 


NGC7624 


61.1 


0.80±0.05 


2.04±0.06 


NGC7625 


24.4 


1.18±0.06 


2.78±0.08 


NGC7673 


48.9 


5.96±0.28 


7.70±0.21 


NGC7678 


49.7 


4.87±0.23 


7.25±0.20 


NGC7798 


34.6 


3.65±0.17 


5.59±0.16 



Table 2. 



Model description 8. 


\ogC F uv h 


logC NU v h 


iogC„ a h 


Cfuv /Chc 


Cnuv /Chc 


O-FUV 


cnuv 


K98 


-43.147 


-42.975 


-41.102 


0.0090 


0.0134 


0.38 


0.21 


Salpeter IMF, Z , 100 Myr 


-43.170 


-42.959 


-41.056 


0.0077 


0.0125 


0.48 


0.24 


Kroupa IMF, Z , 100 Myr 


-43.350 


-43.137 


-41.257 


0.0081 


0.0132 


0.45 


0.22 


Salpeter IMF, Z , 1 Gyr 


-43.207 


-43.023 


-41.056 


0.0071 


0.0108 


0.53 


0.30 


Kroupa IMF, Z , 1 Gyr 


-43.384 


-43.196 


-41.257 


0.0075 


0.0115 


0.50 


0.27 



a For all cases, a constant star formation history and a mass range of 0.1-100Mq are assumed. 

b CFUV> C-nuv and Cu a relate to their corresponding SFRs by the relation SFR(A) = C\ ■ L(A) and are in units of 
Mq yr -1 /ergs s _1 . 



Table 3. Summary of Coefficients 



Relation 


Cocfficient b 


Dispersion 


L(FUV) obs + a*L(TIR) 


0.46±0.12 


0.09 


L(FUV) obs + a*L(25 fan) 


3.89±0.15 


0.13 


L(FUV) obs + a*L 1 . 4G Hz 


71.81±5.06 x 1.E20 


0.14 


L(NUV) obs + a*L(TIR) 


0.27±0.02 c 


0.10 


L(NUV) obs + a*L(25 fim) 


2.26±0.09 c 


0.13 


L(NUV) obs + a*Li. 4G Hz 


41.75±2.97 c x 1.E20 


0.14 



b The coefficients a in the combinations of UV and TIR or 25 fim 
luminosities arc unitlcss, and the luminosities in these combinations 
are in units of ergs s —1 . The coefficients a in the combinations of UV 
and 1.4 GHz radio luminosities are in units of 10~ 7 Hz , the UV lu- 
minosities are in units of ergs s -1 and the 1.4 GHz radio luminosities 
are in units of w Hz^ 1 . 



c The error does not include the uncertainties in the reference lu- 
minosities, so it is under-estimated. 



